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^ ■ The group classification of a class of variable coefficient reaction-diffusion equations with 

I exponential nonlinearities is carried out up to both the equivalence generated by the corre- 

sponding generalized equivalence group and the general point equivalence. The set of admis- 
sible transformations of this class is exhaustively described via finding the complete family 
of maximal normalized subclasses and associated conditional equivalence groups. Limit 
processes between variable coefficient reaction-diffusion equations with power nonlinearities 
■ and those with exponential nonlinearities are simultaneously studied with limit processes 

between objects related to these equations (including Lie symmetries, exact solutions and 
conservation laws) . 



^ . 1 Introduction 
00. 

, The problem of extended group symmetry analysis of variable coefficient reaction-diffusion equa- 

tions of the general form 

in; 

f{x)ut = ig{x)A{u)u,)^ + hix)Biu), (1) 



where fgA ^ 0, was initiated in |23y24j. The case of A and B being power functions, i.e., the 
class of equations having the form 

>< : f{x)ut = {g{x)u''u,), + h{x)u^, (2) 



where fg ^ and {n,mh) ^ (0,0), was successfully investigated. For many reasons, the 
natural continuation of that study is to consider equations of the form ([T]) with A and B being 
exponential functions, 

f{x)ut = (5(x)e™^/,),. + h{x)e"'^. (3) 

Here / = f{x), g = g{x) and h = h[x) are arbitrary smooth functions of the variable x, fg ^ 
and n and m are arbitrary constants. The linear case, which is singled out by the condition 
n = m = 0, is excluded from consideration as it is well investigated. The semilinear equations 
of the form ([3]), which correspond to the constraints n = and m 7^ 0, were already considered 
in [22.i25j . Moreover, equations of the form ([3]) with n 7^ are not related to linear and semilinear 
equations of the same form via point transformations. This is why in the present paper we study 
only the class of equations of the form ([3|) with fgn 7^ 0, which we briefly call class ([3]). Note 
that the parameter n can be gauged to 1 by a simple scaling of variables from the very beginning 
but we will not use this gauge in the present paper and will deal with the general form ([3]). 
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Motivated by the work in [23] as a general outline for similar studies, we carry out the 
complete group classification of class In order to achieve this in the easiest way, we apply 
a number of modern tools of group analysis of differential equations: generalized extended 
equivalence groups, conditional equivalence groups, maximal normalized subclasses, the method 
of furcate split, variable gauges of arbitrary elements by equivalence transformations, additional 
equivalence transformations, etc. [Tn l [T5 l [n i [2n l [23 l [24] . 

The structure of this paper is as follows. At first, in Section [2] using the direct method we 
derive the determining equations for admissible point transformations in class ([3]) . In Section [3] 
we find equivalence groups of different kinds for class ([3]) and show that the consideration of 
this class can be simplified by setting the gauge g = 1- Moreover, two possibilities for further 
simplification in the case m = n via using the associated conditional equivalence group of class ([3]) 
are discussed. In Section S] the group classification of class ([3]) is carried out up to equivalence 
generated by the generalized extended equivalence group of this class. The usage of the method 
of furcate split for solving the group classification problem is explained in detail. Admissible 
point transformations of equations from class ([3]) are exhaustively described in Section[5]in terms 
of conditional equivalence groups and normalized subclasses. Section [6] is devoted to the study 
of contractions (nontrivial limit processes) between equations from classes ([2]) and ([3]). Using the 
derived contractions and the results obtained for class ^ in [23], we construct exact solutions 
and conservation laws for equations from class ([3|). The results of the paper are summed up in 
the conclusion. The thorough gauging of constant parameters in the group classification list for 
class ([2|), which was not presented in [23], is discussed in the appendix. 




2 Preliminary study of admissible transformations 

Due to a special structure of equations from class ([3]), the following problem can be solved 
completely. To describe all point transformations each of which connects a pair of equations from 
class ([3]). Such transformations are called form-preserving [12] or admissible [H] or allowed [26] 
transformations. See [18] for stronger definitions. They can be naturally interpreted in terms of 
the category theory [21]. Note that there exists an infinitesimal equivalent of this notion [2]. 

At first we make a preliminary study of admissible transformations for class ^ using the 
direct method [12]. We briefly describe the corresponding procedure. In what follows we refer 
to ([3]) as a single equation under assuming that the arbitrary elements /, g, h, n and m are fixed 
and as a class of equations if the arbitrary elements are varied. We use the same agreement for 
other classes of equations. Consider a pair of equations from the class under consideration, i.e., 
equation ^ and the equation 



and assume that these equations are connected via a point transformation T of the general form 

i=T{t,x,u), x = X{t,x,u), u = U{t,x,u), 

where \d(T, X,U)/d{t, x,u)\ ^ 0. We have to derive the determining equations for the func- 
tions T, X and U and then to solve them. Simultaneously we have to find the connection 
between arbitrary elements of equations ([3|) and ([4|). 

After substituting the expressions for the new variables (with tildes) into ([!]), we obtain 
an equation in the old variables (without tildes). It should be an identity on the manifold C 
determined by ^ in the second-order jet space with the independent variables (i, x) and the 
dependent variable u. To involve the constraint between variables of on the manifold £, we 
substitute the expression of ut implied by equation ([3]) . The splitting of this identity with respect 
to the derivatives utx, utt, Uxx and Ux implies the determining equations for the functions T, 




/(5)% = {g{x)e''''ux)x + h{x)e 



mu 



(4) 
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X and U . In particular, we obtain that = = = (and hence TtX^Uu 7^ 0). This 
well agrees with results on more general classes of evolution equations [12l[T3l[T?l l21]. The 
restrictions = = are common for point transformations between (1 + l)-dimensional 
evolution equations of order greater than one. The restriction X^ = follows from the additional 
condition that the related equations are quasilinear, i.e., they have the form ut = F{t, x, u)uxx + 
G{t, X, n, Ux) with F ^ 0. 

We take into account the above determining equations. Then collecting the coefficients of Uxx 
gives fgX^e"''^ = fgTte^^. Dividing the last equation by e"" and differentiating the resulting 
equation with respect to u, we obtain that f gTte"'^ {hUu — n) = 0, i.e., hUu = n. Therefore, 
Uu is a nonvanishing constant, which we denote by ^3, and the component U can be represented 
in the form U = 63U + 'ip{t, x) with the smooth function ^ of t and x. The collecting coefficients 
of Ux leads to the equation 

(^2~gVxXx - ^ ^gx + mVX^ - ~g VXxx^ Tte''^ + /X'^t = 0, 

— lb 

where by V we denote c's , which is a nonvanishing function of t and x. The further splitting 
with respect to u implies, in particular, that Xt = 0. Taking into account the found constraints, 
we deduce the following expressions for components of the point transformation T: 

T = T{t), X = ip{x), U = 63U + ^P{t,x), h = -, (5) 

CIS 

where T, ip and ip are smooth functions of their arguments, ^3 is a constant and TfcpxSs 7^ 0. The 
remainder of determining equations for T is given by the unused constraints associated with the 
coefficients Uxx and Ux and the collection of the term without derivatives of n, 

f~gVTt - }g^l = 0, 2^ - ^ + ^4h=x - ^ = 0, (6) 

y g g fx 

— (9—) e"" - - n^/ie™" + ^Tte'^^he'^^^'' = 0. (7) 
(fx \ fxj X y f ^3 

(We additionally simplify the above system by substituting the expression for the ratio /// from 
the first equation into the second one.) 



3 Equivalence groups 

In order to solve a group classification problem it is essential to derive the equivalence transfor- 
mations which preserve the differential structure of the class of differential equations under con- 
sideration and transform only its arbitrary elements. These transformations form a group |17j . 
which is called the equivalence group of the class. 

There exist several kinds of equivalence groups. The usual equivalence group consists of the 
nondegenerate point transformations of independent and dependent variables as well as arbitrary 
elements of the class. Moreover, transformations of independent and dependent variables do not 
depend on arbitrary elements. If such dependence arises then the corresponding equivalence 
group is called generalized [T3]. If new arbitrary elements are expressed via old ones in some 
non-point, possibly nonlocal, way (e.g. new arbitrary elements are determined via integrals 
of old ones) then the equivalence transformations are called extended. The first examples of 
a generalized equivalence group and of an extended equivalence group were presented in [H] 
and [9], respectively. See a number of examples on different equivalence groups and their role in 
solving complicated group classification problems, e.g., in p ^[23t [M]. 
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It appears that the generahzed extended equivalence group G~ of class ([3]) is nontrivial, i.e., 
it is wider than the usual equivalence group of the same class. We construct the group using 
the constraints derived in the course of the preliminary study of admissible transformations 
and varying arbitrary elements in equations Q and ([7|), which form the complete system of 
classifying conditions for admissible transformations in class ([3]). 

For the general values of n and m the exponents e^", e"" and e*"" are linearly independent. 
This is why due to varying of arbitrary elements we necessarily have that rh = m/S^. Splitting of 
equation ([7|) with respect to u results in Vf = (i.e., ip{t, x) = ip{x)) and therefore differentiating 
both sides of the first equation of ([6]) with respect to t we have Tu = 0, i.e., T = 5it + 82- Then 
solving equations ([6]) we obtain 

- _ Spfx 7 _ ^pTt , 

where Sq is a nonvanishing constant. The remaining equations are 

~g^) =0, hfe?^'^5i-5jh = 0. 
Solving these equations, we obtain 

h = ^e-n^K V=(6j'^ + 5X\ (8) 



Vipx ' V J 9{x) 

where (^4 and 6^ are arbitrary constants with {84,, 5^) 7^ (0,0). 

The equivalence group of the subclass singled out from class ^ by the constraint m = n is 
wider than the equivalence group of the whole class ([3]). In other words, the constraint m = n is 
associated with a nontrivial conditional equivalence group of class ([3]). In this case we have less 
splitting with respect io u. As a result, V (or, equivalently, V) is an arbitrary smooth function 
of t and 



nVipx 



V \ y2 



The above results are summarized in Theorems [T] and [2l In what follows we use the notation 

——tb 

^ = 1/V = 6*3 in order to simplify formulas. 

Theorem 1. The generalized extended equivalence group of class ^ is formed by the 
transformations 

i = 6it + 62, x = (p{x), u = 5su + ip{x), 

f = g = 5o<fx-^^g: h = e 's^-^h, n=—, m = — , 

Vx Vx 03 ^3 

where ip is an arbitrary smooth function of x, (fx / 0; -0 and are determined by the formulas 

ij{x) = -^-^\n\^{x)\, ^{x) = 54 Z^^ + '^s. 
n J g{x) 

(5j, j = 0, . . . , 5, are arbitrary constants, 5q5i5z 7^ and (^4, ^5) 7^ (0, 0). 

Thus, elements of G~ are parameterized by six arbitrary constants and a single arbitrary 
smooth function of x. The usual equivalence group G^ of class ^ is the subgroup of the 
generalized extended equivalence group which is singled out with the condition (^4 = 0. 
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Theorem 2. The generalized extended equivalence group of the class of equations 



/(xK = (5(x)e"%,), + Mx)e™ with nfg^O, (9) 
coincides with the usual equivalence group G~^^ of this class and consists of the transformations 
t = dit + d2, x = ip{x), u = 5su + ip{x), 

f = ¥l-^f^ ~g = So^,-i'^g, h=^-^[nh^ + {g^,U^, n = ^, 

where 6j, j = 0,1,2,3, are arbitrary constants, Sq6i6s ^ 0, and -0 are arbitrary smooth 
functions of x with ipx 7^ 0, ^'(x) = e ''3 . 

Elements of are parameterized by four arbitrary constants and two arbitrary smooth 

functions of the single variable x. Therefore, the group is really a nontrivial conditional 

equivalence group of class ([3]). 

In view of Theorem [H the family of equivalence transformations 

— -, u = u, (10) 
xo 9[y) 

parameterized by the arbitrary element g maps class ([3]) onto its subclass consisting of equations 
of the form 

with 5 = 1. The new arbitrary elements are expressed via the old ones in the following way: 

/ = fg, h = gh, m = m, h = n. 

Hence, within the framework of symmetry analysis it suffices, without loss of generality, to 
investigate the equations of the general form 

/(x)txt = (e'^^n,.). + /i(x)e'"" with nf ^ (11) 

instead of class ([3]) because all results on symmetries, solutions and conservation laws of equa- 
tions from subclass (jlip can be extended to the entire class ([3]) with transformations (|10p . In 
other words, up to G"^ -equivalence we can assign the gauge g = 1 for the arbitrary element 50 
The description of the generalized equivalence group of class (jlip and the generalized con- 
ditional equivalence group of the same class which is associated with the condition m = n is 
deduced from Theorems 1 and 2 by setting g = g = 1. 

Theorem 3. The generalized equivalence group Gi of class (|lip is formed by the transforma- 
tions 

t = 0it + 02-, X=- — , U = OsU ln|04X + 05|, 

OiX + 05 n 
f = ^\(>ix + b^\-' f , /i = — |^4x + ^sl n^'^/i, '^ = Y^' 

where 6j, j = 1, . . . ,7 , are arbitrary constants such that 6163 ^ 0, A = 656^ — d^dT ^ and the 
tuple {64,, 65, 65,67) is defined up to a nonzero multiplier; e.g., we can set A = ±1. 



^Different simple gauges of arbitrary elements by equivalence transformations are possible. For example, 
instead of jr = 1 we can set f — 1 but the gauge (? = 1 is more convenient because it results in simpler group 
classification of class ([3]). Simultaneously, we can assign the value 1 to the constant arbitrary element n. It has 
been noted in the introduction that this gauge is not essential in the course of group classification and hence it 
will be used only in the presentation of the final classification list. 
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Theorem 4. The class of equations 



f{x)ut = (e""u^)^ + /i(x)e"" with nf ^ (12) 
admits the generalized equivalence group G'J"^^^ consisting of the transformations 

i=5it + 52, x = ip{x), u = 63U+ —ln\6l(px\, (13) 

2n 

7 J J- I 1-3 „ I „ _2, , 1-3, 1 _ n 

n 

where 6j, j = 0, ... ,3, are arbitrary constants with S06163 7^ and (/? = (p{x) is an arbitrary 
smooth function with ipx ^ 0- 



Class p2p can be mapped onto a proper subclass with only one arbitrary element depending 
on X using an appropriate family of point transformations from the group G^^^^. The most 
convenient gauges for arbitrary elements are the gauges / = 1 and h = 0. 

The first gauge can be realized by the transformation 

r 2 1 

t = sign{f)t, x= f{y)^dy, u = u + —\n\f\, (14) 

J Xq 

which maps an equation of the form (jl2p to the equation ui = {e'^^Ux)x + ^e"", i-e., / = 1- The 
other new arbitrary elements are expressed via the old ones in the following way: 

h = j(^f'^h + ^{f-^i)x^, h = n. 

Theorem 5. The generalized equivalence group Gj^g^i of the class of equations 

ut = (e^^-u^)^ + /i(x)e"" with n / 0, (15) 
consists of the transformations 

i=6it + 62, x = SiX + 65, ■u = 53M + — In^, h = ^h, n = 

n di 01 03 

where 6j, j = 1, . . . ,5 are arbitrary constants, did^di ^ 0, 5i > 0. 



The gauge h = for class (jl2p can be realized by the family of transformations ()13p . where 
Si = 5^ = 1, 82 = and the function (p satisfies the ODE (Iv'zl"^ )2:2: + nh\ipx\~^ = 0. 

Theorem 6. The generalized equivalence group G^f^^Q of the class of equations 

/(a;)ut = (e""u^)^ with n/ / 0, (16) 

is projection of the group Gi on the space {t, x,u, f,n). 

It is possible to carry out the complete group classification of equations from class (|12p using 
either the gauge / = 1 or the gauge h = 0. The choice of the first gauge is justified in Section HI 
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4 Lie symmetries 



In general, the group classification of class ([3]) is carried out within the framework of the classical 
Lie approach |16tll7|. At the same time, we additionally apply a number of modern tools of 
symmetry analysis, which essentially simplify both related calculations and the presentation of 
final results. Thus, gauging the arbitrary element 5 to 1 by equivalence transformations, we in 
fact classify subclass (jlip instead of the entire class. The main equivalence relation involved in 
the consideration is generated by the generalized equivalence group of this subclass, which 
contains the usual equivalence group of the same subclass as a proper subgroup. In other words, 
we use G^^-equivalence, which is stronger than -equivalence prescribed by the classical Lie 
approach. Moreover, for group classification of equations from subclass (jlip with m = n we 
involve the equivalence relation which is generated by the conditional generalized equivalence 
group Gij^^^ and is even stronger than G^" -equivalence. In total, this leads to the reduction 
of classification cases and lowering the number of additional equivalence transformations to be 
constructed. Lie symmetry extensions are separated using the method of furcate split [lOpiSj. 

The reduction of group classification of class ([3]) to that of its subclass ([TT]) is possible due to 
combining two facts. Namely, 1) each equation from class ^ is mapped by a transformation from 
the generalized extended equivalence group of this class to an equation from subclass (fTTI) 
and 2) the generalized equivalence group Gi of class ()lip only consists of those transformations 
from G~ which preserve the constraint g = 1. The same remark is true for class ([9|) with the 
usual equivalence group G~=„ and class (fT2]) with the generalized equivalence group GJ^^^^. 

We look for vector fields of the form 

Q = T{t, X, u)dt + i{t, X, u)dx + r/(t, x, u)du 

which generate one-parameter groups of point symmetry transformations of an equation C from 
class (jlip . These vector fields form the maximal Lie invariance algebra A^^^ = A^^^{C) of the 
equation C. Any such vector field Q satisfies the criterion of infinitesimal invariance, i.e., the 
action of the second prolongation Q^^^ of Q on equation (jlip results in the conditions being an 
identity for all solutions of this equation. Namely, we require that 



identically, modulo equation pT]) . 

After the elimination of ut due to (|lip . equation (jl7p becomes an identity in six variables, 
t, X, u, Ux, Uxx and utx- In fact, equation p7p is a multivariable polynomial in the variables 
Ux, Uxx and utx- The coefficients of different powers of these variables must be zero, giving the 
determining equations on the coefficients r, ^ and rj. Solving these equations, we immediately 
find that r = T(t), ^ = S,{t,x). This completely agrees with the general results on point trans- 
formations between evolution equations [12j. Then the remaining determining equations take 
the form 



Splitting of the first equation with respect to u and the subsequent integration imply that 




(17) 



6/ 



+ ({/Ji + (2t + (m - n)r,) h) e"" 




where 



■rf = rp{i) is a smooth function of t. 
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Finally we obtain the classifying equations which essentially include both the residuary un- 
certainties in the coefficients of the vector field Q and the arbitrary elements of the class under 
consideration: 

?^ = Tt-^C. + nr?o, (18) 

In order to find the common part of Lie symmetries for all equations from class pip (resp. 
class (l3|)), we should vary the arbitrary elements and split with respect to them in the determining 
equations. This results in rt = ^ = = 0. 

Proposition 1. The kernel algebra, i.e., the intersection of the maximal Lie invariance algebras 
of equations from class (|lip (resp. class ^) is the one- dimensional algebra A'^ = {dt). 

The classification of possible extensions of A'^ is reduced to the simultaneous solution of 
equations (|18p and (|19p with respect to both the residuary uncertainties in the coefficients of 
the vector field Q and the arbitrary elements up to G^^-equivalence. Further splitting with 
respect to u in the last equation depends on the value of m. Namely, we should separately 
consider three exclusive cases 

1) rn-T^O, n; 2) m = 0; 3) m = n. 

The case /i = is special as the value of m is undefined in this case but it can be included to 
the case m = n (see Section 2). This is why in the first two cases we assume that h ^ 0. 

4.1 General case of m 

In the first (general) case splitting of (fT9]) gives the equations r/? = and S^xxx = 0. Differentiation 
of (fT8|) with respect to t leads to the equation tu = 0. Therefore, 



T = C4^t + C5, ^ = nC2x'^ + CiX + Co, T] = C2X + C3. (20) 

The classifying equations take the form 



fx 

{nc2X^ + cix + co)-j- = -3nc2X + nca + C4 - 2ci, (21) 
{nc2x'^ + c\x + co)-^ = — (3n + m)c2X + (n — m)c3 — 2ci. (22) 



Further investigation will be carried out using the method of furcate split |10pi5l. For any 
operator Q from A"^^^ the substitution of its coefficients into equations (j2ip and (j22p gives some 
equations on / and h of the general form 

[nax^ + hx + c)^ = — 3nax + d, (23) 

(nax"^ + bx + c)— = — (3n + m)ax + (n — m)p — 2b, (24) 
h 

where a, b, c, d and p are constants which are defined up to nonzero multiplier. The set V of 
values of the coefficient tuple (a, b, c, d, p) obtained by varying of an operator from A"^^^ is a 
linear space. Note that the first three components of any nonzero element from V form a nonzero 
triple (otherwise the corresponding system with respect to / and h would be inconsistent). 
The dimension k = ^(^4™'^^) of the space V is not greater than 2 otherwise the corresponding 



8 



equations form an incompatible system with respect to / and h. The value of k is an invariant 
of the transformations from G^. Therefore, there exist three G^^-inequivalent cases for the value 
of k: k = 0, k = 1 and k = 2. We consider these possibilities separately. 

I. The condition k = means that (j^T|) - ([2^ are not equations with respect to / and h but 
identities. Therefore, / and h are not constrained and Cj = 0, z = 0, . . . , 4. We obtain nothing 
but the kernel algebra A'~^ presented by Case 1 of Table 1. 

II. If A; = 1 then we have, up to a nonzero multiplier, exactly one system of the form (|23p - (|24p 
with respect to the functions / and h. Integrating this system up to G J" -equivalence gives 

f f —3nax + d \ f I (3?t- + m)ax + 26 + (m — n)p \ 

/ = exp / 2-TT — — , /i = eexp - / ^r—- — ■ dx , 

\ J nax'^ + ox + c J \ J nax'^ + ox + c J 

where e = ±1 mod G'^ , and the coefficients of Q take the form 

T = (d + 26 — np)\t + C5, ^ = {nax^ + bx + c)A, r] = {ax + p)X. 

Here A and C5 are arbitrary constants. Hence A^^^ is the two-dimensional algebra given by 
Case 2 of Table 1. 

It is necessary to investigate equivalence of cases with the arbitrary elements / and h of the 
above form, corresponding to different values of the tuple {a,b,c,d,p). 

Lemma 1. Up to G"^ -equivalence the parameter tuple {a,b,c,d,p) can be assumed to belong to 
the set {{0,1,0, d,p), (0,0,l,l,p), (0,0,1,0,1), (0,0,1,0,0), (1/n, 0, 1, d,p)}, where d^ -3/2 
and, if d = —3/2, np ^ 1/2; d ^ and, if d = 0, p 0. 

Proof. Combined with the multiplication by a nonzero constant, each transformation from the 
equivalence group G^ is extended to the coefficient tuple of system (j23|) -- (pl]) : 

hd = v{5lna — S^S^b + d^c), 6 = v{—25^5'jna + {SiS'j + 5^5Q)b — 254,5qc), 
c = v{djna — 6e6jb + 6qc), d = uAd + 3i'{656jna — d^djb + S^dQc), 
hp = uAnp — ^{d^djna — S^djb + 64,6qc), 

where A = 6^6^ — 54^6j 7^ 0, is an arbitrary nonzero constant. 

There are only three Gj'-inequivalent values of the triple (a, b, c) depending on the sign of 
D = b'^ — Anac, 

(0,1,0) if D>0, (0,0,1) if D = 0, (l/n,0,l) if D<0. 

Indeed, if > then there exist two linearly independent pairs {84,65) and {6q,6j) such that 
S^na — S^S^b + (5|c = and 6jna — SeSjb + d^c = 0. For these values of (5's we have d = c = 0. 
In the case D = we choose values of 64, 65, 5q and 67 for which — 6465b + 6lc = 
and the pair {5g,5j) is not proportional to the pair (^4,^5). Then we obtain that a = and 
b = u67{64b — 265na) + z/56(<^5^ — 264c) = 0. The residual coefficient (6 if D > and c if D = 0) 
is necessarily nonzero and hence can be scaled to 1 by choosing the appropriate value of i^. If 
D < 0, we have ac^ and can set a > 0. The matrix 

/ na -6/2 \ 

v-6/2 c ; 

is symmetric and positive. Hence the corresponding bilinear form is a well-defined scalar product. 
We choose u = 1 and pairs {64,65) and {6^,67) which are orthonormal with respect to this 
product. Then hd = c = 1 and 6 = 0. 
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Certain freedom in varying group parameters is preserved even after fixing one of the above 
inequivalent forms for both the tuples (a, b, c) and (a, b, c). This allows us to set constraints for 
the coefficients d and p. 

Thus, it follows from the equality {a,b,c) = {a,b,c) = (0,1,0) that 6465 = 6q6j = and 
6/i6j + 656Q = . There are two cases for the solution of the above system in 5^s: either 
(54 = ^7 = and A = 5^5% = or (^5 = = and A = —S^Sj = —u^^. In the first case the 
coefficients d and np are identically transformed. In the second case the transformation takes 
the form d = —d — 3, hp = —np + 1. This is why up to G^^-equivalence we can assume that 
d ^ -3/2 and, if d = -3/2, hp ^ 1/2. 

The equality (a, b, c) = (a, b, c) = (0, 0, 1) implies ^4 = 0, i^5q = 1 and A = 656^ / 0. 
The transformation of d and np is reduced to simultaneous scaling with the same multipliers, 
d = uAd and hp = lyAnp. This allows us either to set d = 1 if d 7^ or to scale p if d = 
and hence d = 0. Therefore, we obtain the following inequivalent tuples (0,0, 1, l,p), if d 7^ 0, 
(0, 0, 1, 0,1) if d = 0, p / 0, and (0, 0, 1, 0, 0) in the case d = p = 0. The last tuple corresponds 
to the constant-coefficient equations (jlll) (Case 3 of Table 1). 

Setting [na, b, c) = (no, b, c) = (1,0, 1) results in 6^ + 5'^ = + 6j = u^^ and 646Q + d^dj = 0. 
Hence 6q = eS^ and 5-/ = —£64, where e = ±1. The transformation of the coefficients d and np 
is reduced to the multiplication by e, d = ed and hp = enp. This is why we can only set d ^ 
and, if d = 0, p ^ 0. □ 

Lemma [T] implies that up to G^^-equivalence the case A; = 1 is partitioned into the three 
inequivalent subcases: 

1. (/,/i) = (|x|^e|x|9), q = {n-m)p-2: {du {d + 2 - pn)tdt + xd^ + pdu); 

2. (/, /i) = (e'^^,ee''^), q = {n - m)p: {dt, {d - pn) tdt + + pdu); 

3. (/,/i) = (^(x2 + 1)-Ie'='^'''=*'^'^^,e(x2 + l)-|-5Se«^''=*^'^^) , q={n-m)p: 
{dt, {d - pn)tdt + {x^ + l)5x. + {x/n + p)du). 

Here (d, g) 7^ (0,0),(— 3, — 3 — m/n) and (d, g) 7^ (0,0) for the first and second subcases, re- 
spectively. (Otherwise k = 2, see below.) It additionally follows from Lemma [T] that up to 
GJ" -equivalence we can set certain constraints for the parameters d and q. (It is convenient 
to use q instead of p as a parameter.) For the first subcase an exhaustive gauge implied by 
G J" -equivalence consists of the inequalities d ^ —3/2 and, if d = —3/2, q ^ —3/2 — m/{2n). 
They can be set using the equivalence transformation 

1 1 , , , 

t = t, X = —, u = u m X , (25) 

X n 

whose extension to the parameters d and q is given by d = — d — 3 and q = —q — 3 — m/n. In 
the second subcase the parameters d and q can be gauged using a scaling of x. More precisely, 
d = 1 mod Gi if d 7^ and q = 1 mod if d = 0. In the last subcase we can just simultaneously 
alternate the signs of d and q. Hence the exhaustive gauge is presented by d ^ and, if d = 0, 

III. Let k = 2. We choose a basis {(a^, fej, Cj, di,pi), i = 1, 2} of the space V of tuples (a, b, c, d,p) 
associated with and introduce the notation 

Pi = naiX^ + biX + Ci, Qi = —3naiX + dj, Ri = — (3n + m)aiX + (n — m)pi — 2bi. 

Then Pj 7^ 0, i = 1, 2, and the system for the functions / and h is written in the form 

fx _ Qi^ _ Q2_ hx_ _ Ri_ _ R2 

J-p[-^' T-p[-W ^^^^ 
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Consider two different cases depending on whether the first components of all elements from 
the space V are zero. 

If this is true, we can choose a basis of V with bi = C2 = 1 and 62 = ci = 0. Then system (j26|) 
obviously implies that fx = hx = 0. As fh 7^ 0, up to G^" -equivalence we obtain Case 3 of 
Table 1. 

If the space V contains tuples with nonzero first components, we can assume without loss 
of generality that ai = 1 and 02 = 0. Then we have 62 7^ as otherwise system (j26p is not 
compatible. We set 62 = 1 and 61 = by changing the basis of V and then set C2 = using a 
shift of X allowed by Gj'-equivalence. Equations (j26p are consistent if and only if ci = 0. They 
imply that (/, /i) = {x~^ , e\x\^'^^ ~^ ) mod GJ", where e = ±1. The corresponding maximal Lie 
invariance algebra 

{dt, nx'^dx + xdu, {m + n)tdt + (m — n)xdx — 2dy) 

is a two-dimensional extension of the kernel algebra A'^ = {dt). This Lie symmetry extension 
is not included in Table 1 because it is reduced to Case 3 by the transformation (p5|) from the 
equivalence group Gi- 

4.2 Case m = 

If 771 = then ^ = nc2X^ + cix + cq and the classifying equations have the form 

{nc2x'^ + cix + co)^ = n - 3nc2X - |ci -|- nrj^, 

f 3 \ 

Vtf = {nc2X^ + cix + co)hx + i 3nc2X + -ci - nrp \ h. 

Differentiating the latter equation with respect to t, we obtain the equation frj^^ = —nhr]^. If 
{h/f )x 7^ then 77° is a constant and hence the corresponding equation possesses the same 
maximal Lie invariance algebra as in the case of general value of m for the same values of / 
and h. In what follows we assume that {h/f)x = 0, i.e., e = h/f is a constant which equals ±1 
modulo Gi . Substituting the condition h = ef into (j27p . we obtain that 

T = C4e~'^^^ + C5, ^ = nc2X^ + cix + cq, rj = C2X + ec^e''^'^^ + C3, 

and the classifying equation on the function / is the following one 

fx 

(nc2X^ -|- cix + co)-j- = —3nc2X — 2ci -|- nca. 

It is obvious that this equation has the form (j23p . Applying the method of furcate split in the 
same way as in Section 14. H we obtain only three G^'-inequivalent Lie symmetry extensions, 
namely. Cases 4-6 of Table 1. Case 5 is partitioned into the three G^^-inequivalent subcases: 

1. f = \x\'i: {dt,e-'''\dt + edu),nxdx + {d + 2)du); 

2. f = : {dt, e-'^\dt + ed^), ndx + 

3. / = + i)-|e'^'^'-cta°^: {dt, e-^"*(9t + edu), n{x^ + l)dx + {x + d)du). 



Here always h = ef. In the first subcase we should assume that d 7^ —3,0 as otherwise this 
subcase is reduced to Case 6 with a wider Lie invariance algebra, and d ^ —3/2 mod G'^, 
cf. Lemma [T] and the transformation (|25p . In the last subcase we can just alternate the sign 
of d. Hence the exhaustive gauge up to G^" -equivalence is presented hy d ^ 0. 
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4.3 Case m = n 



If m = n then rj is a constant and the classifying equations take the form 

ih = n- ^ix + nrf, ^i^^^ + + 2Cxh = 0. (28) 
J 2 In 

To carry out the group classification in this case, it is necessary to use an additional gauge of 
the arbitrary elements of class (fT2|) . Consider both the gauges proposed in Section 2, namely, 
/ = 1 and /i = 0. 

Choosing the gauge / = 1 results in subclass p^ . The first equation of (f28|) implies that 
ixx = 0, i.e., ^ = cix + C2, where ci and C2 are constants and r] = -^ix + ri^ \s a, constant, which 
is denoted by C3. Then r = (2ci — ncs)^ + C2. The classifying equation on the function h can be 
written in the form 

(cix + C2)hx + 2cih = 0. 

It is easy to derive Cases 7, 8, 9 and 10 of Table 1. 

Introducing the gauge /i = 0, we replace the study of class (fT2]) by the study of class (fT6]l . 
Under this gauge, determining equations gives ^xxx = and hence r, ^ and rj have the form ([20]) . 
The classifying equation on the function / coincides with equation (12ip . The solution of (I2ip 
results in the following cases of Lie symmetries extensions: 

arbitrary /: ntdt - (9„); 

/ = exp I I 3naa; + d ^\ ^ ^^^^ ^^^^ _ 26)t9t + (nax^ + 6x + c)5a; + ax5„); 

\J nax'^ + ox + c J 

/ = 1: {dt, ntdt - du, d^, 2tdt + xd^). 

The first and third cases obviously correspond to Cases 7 and 10 of Table 1, respectively. 

Similarly to the general case of m and the case m = 0, the second case of Lie symmetry 
extensions with m = n under the gauge h = is partitioned into the following three Gi- 
inequivalent subcases: 

1. f = \x\^: {dt,ntdt-du, {d + 2)tdt + xd^); 

2. / = : {dt, ntdt - <9„, tdt + d^); 

3. / = + i)-|e'^'^'-ctan^' : {dt, ntdt - du, ndtdt + n{x^ + 1)9^. + xdu). 

In the first subcase we should again assume that d 7^ —3, as the value d = corresponds to 
the case of / = 1 and h = with a wider Lie invariance algebra (Case 10 of Table 1) and the 
value d = —3 is G^" -equivalent to the value d = 0, cf. the consideration for A; = 2 in Section HTTl 
Additionally, Lemma [J implies that using the transformation ()25p we can set the gauge and 
d ^ —3/2 mod G^. In the last subcase the equivalence transformation alternating the sign of x 
leads to the gauge d ^ mod Gi . 

Remark 1. The above subcases are related to cases of Table 1 via point transformations of 
the form (flil) which belong to the equivalence group G^^^^. Thus, the first subcase with 
d {—3, —3/2,0}, the second subcase and the third subcase with d ^ are mapped to Case 8 
with 

d{d + 3) 1 , d'^ + 9 

a = , , „ , a = — and a = — — ^5-, 
4n(d + 3/2)2' 4^ 4^(^2 

respectively. The first subcase with d = —3/2 and the third subcase with d = are mapped to 
Case 9 with e = — l/(4n) and e = 1/n, respectively. 
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Analyzing the group classifications of the class (fT2|) under both the gauges, we conclude 
that each of the gauges / = 1 and h = has certain advantages and disadvantages. More 
precisely, under the gauge h = equations have a less number of summands and the group 
classification list in the case m = n is similar to the group classification lists in the cases of 
general m and m = 0. The gauge /i = is also convenient in order to look for additional 
equivalence transformations between the cases m = and m = n. At the same time, the 
gauging transformation, the corresponding equivalence group and the inequivalent values of the 
residuary arbitrary element, h, arising in the course of group classification are much simpler for 
the gauge / = 1. In particular, it is very easy to partition into inequivalent cases under this 
gauge. Moreover, the group classification list presented in Table 1, which involves the gauge 
/ = 1 in the case m = n, is well consistent with the group classification list found for class ([2]) 
in [23]. This is important for the study of contractions between cases of Lie symmetry extensions. 

4.4 Classification list and additional equivalence transformations 

Results of group classification for subclasses of class ([3]) are collected in Table 1. 



Table 1. Results of group classification of class ((3]) under the gauge g — 1. 



no. 


/(^) 


h{x) 


Basis of yi™'"' 


General case of m 


1 


V 


V 


dt 


2 




hi{x) 


dt, (d + 26 - pn)tdt + {nax^ + bx + c)dx + {ax + p)du 


3 


1 


e 


dt, dx, 2mtdt + (m — n)xdx — 'idu 


m = 0, h^O, {h/f)x = 


4 


V 


ef{x) 


dt, e-''^\dt + edu) 


5 




efi{x) 


dt, e~^"\dt + edu), n{nax^ + hx + c)dx + {nax + 2b + d)du 


6 


1 


e 


dt, e-="*(9t + edu), dx, nxdx + 2du 


m — n or h = 


7 


V 


V 


dt, ntdt ~ du 


8 


1 


ax 


dt, ntdt - du, nxdx + 2du 


9 


1 


e 


dt, ntdt - du, dx 


10 


1 





dt, ntdt - du, dx, nxdx + 29„ 



Here n, a and £ are nonzero constants, n = 1 mod Gi , e = ±1 mod G^, 

and up to G^-equivalence the parameter tuple (a, b, c, d, p) can be assumed to belong to the set 

{(0,l,0,d,(g + 2)/(n-m)), (0,0,l,l,p), (0,0,1,0,1), (1/n, 0, 1, d,p)}, 

where (d, q) (0, 0), (—3, —3 — m/n) and modulo Gi we can also set d ^ —3/2 and, if d = —3/2, q ^ —3/2 — 
m/(2n); d ^ and, if d = 0, p ^ 0. In Case 5 the parameter p should be neglected. In Case 7 the arbitrary 
element / (resp. h) can be additionally gauged by transformations from Gi „^^„. For example, we can set f — 1. 

There exist additional equivalence transformations between classification cases presented in 
Table 1. Namely the point transformation 

t' = — e="*, x' = x, u' = u-et (29) 
en 
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links the equations f{x)ut = {g{x)e^^Ux)x +£f{x) and f{x')u^, = {g{x')e^'^' u'^,)^, . This trans- 
formation belongs to no equivalence group found in Section 2 and reduces Cases 4-6 of Table 1 
to the set of cases 'm = n or /i = 0'. For the reduction to be precise to Cases 7-10 of Table 1, 
for Case 5 transformation (j29p should be composed with an appropriate transformation of the 
form (dH). Such compositions map subcases of Case 5 to subcases of Case 8 and 9, cf. Remark[TJ 
The transformations described exhaust additional equivalence transformations within the clas- 
sification list from Table 1. It is proved in the next section within the framework of admissible 
transformations. Thus, transformation (I29p can also be included in the framework of conditional 
equivalence but the corresponding conditional equivalence group is too complicated. 
As a result, we obtain the following assertion. 

Theorem 7. Up to point transformations, a complete list of Lie symmetry extensions for equa- 
tions from class ([3]) is exhausted by Cases 1-3 and 7-10 of Table 1. 

Therefore, the maximal dimension of Lie symmetry algebras for equations from class ([3]) 
equals four, and for equations with m ^ 0,n this dimension equals three. 

In the Table 1 we do not present the partitions of Cases 2 and 5 into inequivalent subclasses 
in detail. The complete description of these partitions is given in Sections 14.11 and l4.2[ 

Corollary 1. If an equation from class (0) is invariant with respect to a four- dimensional Lie 
algebra then it is reduced using point transformations to the equation Ut — (^e^Ux^x • 

Corollary 2. If an equation from class f^j with m ^ 0,n possesses a three-dimensional Lie 
invariance algebra then it is mapped by a point transformation to the equation ut = {e^Ux)x^e"^^ ■ 



5 Classification of admissible (form-preserving) transformations 

In contrast to the construction of equivalence groups, in the course of the complete study of 
admissible transformations we cannot vary arbitrary elements. The description of the set of 
admissible transformations will be presented in terms of conditional equivalence groups and 
normalized subclasses. 

The consideration of Section [2] implies that any admissible transformation in class ([3]) is of 
the general form ([5]), where the parameters satisfy the inequality Ttfx^s 7^ and equations ^ 
and ([7]). Differentiating the first equation of ([6]) with respect to t, we have that V = 9{x)/Tt. 

— lb 

Recall that V denotes the expression e^3 . Then the determining equations become 

fgG - fg^x = 0, 2— + —Lpx = 0, 30) 

g g <fx 

— [g—] e"" - fTt (—] - n^he"^^ + ^Tte^^he'^^'^'' = 0. (31) 



Vx \ Vx / X \^t/t J "a 

Solving equations (j30p . we express the new arbitrary elements / and g via the old ones: 



f So . ^ Soifx 
J = TT— /' 9 = -0^9, 



where Sq is an arbitrary nonzero constant. 

The splitting of equation ([3T]) with respect to u and the subsequent integration of the deter- 
mining equations appreciably depend on values of m and rh. Consider possible cases for these 
values. 
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If m = m/(53 7^ 0, equation ([HT]) implies Tj^ = 0, i.e., T = 5it + 82 and therefore ^ = 9{x)/5\. 
The other conditions obtained from (j31|) result in 

m = — ^ U,n: n = — — — e *3 ^ /i, where & = I 04 / — — + 05 
03 oidifx V ./ 5(2;) 



m n ~ 6063 
m = — = —: n 



nh ( gOx 



where 6 is an arbitrary function. 



^3 82. ' n6(px 

Therefore, in the case m ^ 0,n (resp. the case m = n) any admissible transformation is induced 
by a transformation from the generalized extended equivalence group G"^ (resp. the conditional 
equivalence group G~^„). 

If m = m = then equation (|3ip implies an expression for 9 and an equation connecting the 
ratios h/ f and h/ f, 

In view of the simplest differential consequence 



of the last equation we have {h/ f )x 7^ and Tu = if {h/ f )x 7^ 0. Therefore, any admissible 
transformation for equations with m = and {h/f)x 7^ is induced by a transformation from 
the equivalence group . 

This is not the case if {h/f)x = 0. Then {h/f)x = 0, i.e., the condition {h/f)x = is invariant 
with respect to Gi- We denote the constants h/f and h/f by a and q, respectively. Then we 
get the following equation for the function T{t): 

l\ 1 

— = —na— + na, 
TtJt Tt 

where n = n/d^. We integrate this equation and present the general solution in such a form that 
continuous dependence of it on the parameters a and a is obvious: 

gHoT _ -j^ ^nat _ e"°"^ — 1 

00 7^0: — = Si 1-62, a = 0, 07^0: — = 5it + 82, 

na na na /onx 

^nat 

07^0, a = 0: T = 5i 1-52; a = a = 0: T = 5it + 82- 

na 

Here 5i 7^ as Tt 7^ 0. 

Theorem 8. The generalized extended equivalence group G^^q ^^^^^ of the subclass of 
class ([3]), which is singled out by the conditions m = and {h/f)x = 0, consists of the transfor- 
mations 

t = T{t), x = ip{x), u = 63U + 7p(t, x), 

f = —^f, g = 6o^x^^g, n = —, 

where the smooth function T = T{t) with Tt ^ is defined by (j32p . a = h/f and a = h/f are 
constants, ip is an arbitrary smooth function of x with (fx 7^ 0, the functions ip and are defined 
by the formulas 

^l^{t,x) = -^-^ln\Ttm{x)\, ^(x) = J4 / 4^ + -^s, 
n J g{x) 

6j, j = 0, . . . ,5, are arbitrary constants, ^o'^i'^s '^^'^ ('^4) ^5) 7^ (0) 0)- 
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In contrast to , we do not use G~ j^^y in the course of group classification of 

class pip because the application of this conditional equivalence group does not have a crucial 
influence on classification, and the corresponding system m = 0, {h/ f)x = for arbitrary 
elements is less obvious. At the same time, transformations from G^^q ^^^^^ play the role 
of additional equivalence transformations after completing the classification (see the previous 
section). 

If /i/i 7^ 0, m = and m = h then the determining equations yield 

Hence both equations ([3]) and @ are mapped to the subclass '/i = 0' of the class under consider- 
ation. These mappings are realized by the transformations from the corresponding conditional 
equivalence groups. We can assume that their images coincide. Therefore, the admissible trans- 
formation is composition of a transformation 71 of ([3|) to the equation f{x)ut = {g{x)e"'^Ux)x 
and a transformation Ti of the equation f{x)ut = {g{x)e'^^Ux)x to The transformation 71 
belongs to G'^^q (h/f) =o transformation Ti belongs to G~^„, and we can choose ip = x 

in 71. 

The case hh ^ 0, m = n and m = is considered in similar way. 

All the possible cases are exhausted. We summarize the investigation of admissible transfor- 
mations in class ^ in the following assertion. 

Theorem 9. Let the equations 

f{x)ut = {g{x)e''''ux)x + /i(x)e™" and j{i)ut = (5(x)e""^/,), + /i(x)e'^" 

he connected via a point transformation T in the variables t, x and u. Then 

either — = — or {m^fh) = (0,n) or (m,fh) = (n,0). 
h n 

The transformation T is induced by a transformation from 

a) if either m ^ 0,n or m = 0, {h/ f)x 7^ 0; 

b) Gm=n if m = n and m / 0, then also rh = h; 

c) (^y^j^^Q if m = fh = 0, {h/f)x = 0, then also {h/f)x = 0. 

If m = and rh = n then {h/f)x = and the transformation T is the composition of two 
transformations, from G^^q ^^^^^^ and with the intermediate equation having h = 0. 

The case with m = n and rh = is similar to the previous one. 

Corollary 3. Class Q is represented as the union of its three maximal normalized subclasses 
separated by the conditions 

{h^O, m / 0, n) or (m = 0, {h/f)x / 0); m = 0, {h/f)x = 0; m = n. 

Only the latter two subclasses have a non-empty intersection, and the intersection being the 
normalized subclass 'h = 0'. 

Recall that the class of differential equations is called normalized if any admissible transfor- 
mation in this class belongs to its equivalence group. See [l8] for strong definitions. 

The subsets of equations appearing in Corollary [3] are really subclasses of class ([3]) since they 
are singled out from class ([3]) by usual systems of differential equations and/or inequalities with 
respect to arbitrary elements. This is not obvious only for the first subclass. The condition 
corresponding to it is in fact equivalent to a single inequality, (m — n)h{m? + {{h/f)x)'^) / 0. 
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6 Contractions 

Examples of nontrivial limits between equations admitting Lie symmetry extensions are known 
for a long time. For instance, in [1] equations with exponential nonlinearities were excluded 
from the group classification list of nonlinear diffusion equations as a separate case and were 
just considered as a limiting case of equations with power nonlinearities. At the same time, it 
looks more convenient to include such cases to classification lists and then indicate connections 
between different classification cases via limiting processes. Using the analogy with theory of 
Lie algebras such connections are called contractions. A theoretical background on contractions 
of differential equations, their Lie symmetry algebras and solutions was first discussed in 

In this section we relate, via contractions, the group classification lists obtained for classes ^ 
and ([2]) in the present paper and in [23j, respectively. For convenience of the presentation, in- 
equivalent cases of Lie symmetry extension in class ([2]) are collected in Table 2. Necessary 
explanations on the equivalence involved and thorough gauges of involved parameters by equiv- 
alence transformations are given in the appendix. Then contractions are used to construct exact 
solutions of equations from class ([3]) using known solutions of equations from class ^ . In Sec- 
tion [63] we demonstrate three different ways of a similar consideration for conservation laws: in 
terms of contractions of associated characteristic or conserved vectors or divergence expressions 
themselves. 



6.1 Contractions of equations and of Lie invariance algebras 

At first we apply the equivalence transformation 

i = 6t, X = Vdx, u = 6{u — 1), h = —, m = — (33) 



parameterized by a positive constant parameter 6 to the equation from class ([2]) with the values 
arbitrary elements g = 1 and / and h presented in Case 2 of Table 2. The constant parameters 
a, b, c, d and p are transformed in the following way 

d = a, b = c = c, d = p = VSp, a = 5a (34) 

wherever this is relevant, i.e., we change parameters if and only if they appear in the values 
of arbitrary elements of the initial equation. Then, we take the imaged equation and proceed 
to the limit 6 — t- +oo. This results in the equation from class (jlip with the values of the 
arbitrary elements / and h presented in Case 2 of Table 1. The same procedure establishes a 
contraction between the associated Lie algebras of vector fields. The corresponding notation will 
be 2.2 — >• 1.2, where the first numbers indicate the numbers of the tables and the second numbers 
indicate the numbers of cases within the tables. We present the complete list of contractions 
which replace power nonlinearities by exponential ones and, therefore, connect cases of Lie 
symmetry extensions for classes ^ and ([2]): 

2.1^1.1, 2.2^1.2, 2.3^1.3, 2.4^1.4, 2.5^1.5, 
2.6^1.6, 2.8^1.7, 2.9^1.8, 2.10^1.9, 2.11^1.10. 

6.2 Contractions of Lie reductions and exact solutions 

In [23j we carried out Lie reductions and constructed Lie exact solution for equations from 
class ^ with the values of arbitrary elements presented in Cases 9 and 12 of Table 2, which 
admit three-dimensional Lie symmetry algebras. It is shown in the previous subsection that 
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Table 2. Results of group classification of class ((2| under the gauge g = 1 [23] 



no. 


n 




h(x) 


Basis of A"""" 


General case 


1 


V 


V 


V 




2 


V 


/l(^) 


/ii(a;) 


9t, (d + 26 — pn)tdt + ((n + l)aa;^ + &a; + c)dx + {ax + p)udu 


3 


V 


1 


e 


dt,dx,2{l — m)tdt + (1 + n - nijxdx + 2udu 


m=l, h^O, {h/f)^=0 


4 


V 


V 


£/ 


dt,e''"'{dt+eudu) 


5 


V 


/l(^) 




at,e"^"'(at +£M5„),n((n + l)ax'^ + bx + c)dx + {nax + 2b + d)udu 


6 


3 


1 


£ 


dt,dx,e~''"\dt + eudu), nxdx + 2u9„ 


7 


4 
3 


1 


£ 


9t, dx, e^^'(9t + eudu), -^xdx + 2udu, —^x'^dx + xudu 


m — n + 1 or h — 


8 


V 


V 


V 


dt , ntdt — udu 


9 


+ 


1 




dt, ntdt - udu, 2tdt + xdx 


10 


+ -\ 


1 


£ 


dt,ntdt ~ udu, dx 


11 


+ -\ 


1 





dt,dx, ntdt - udu , 2tdt + xdx 


12 


4 
3 




Q 


dt,tdt + judu,dx - judu 


13 


4 
3 


1 





dt,dx, ^tdt + udu, 2tdt + xdx, -^x'^dx + xudu 



Here a is arbitrary constant, a 7^ in Case 9, £ = ±1, 
-{3n + 4)ax + 



fi {x) = exp 



(n + l)ax^ + bx + c 



dx 



hi{x) = exp 



-(3(n + 1) + m)ax + {n - m + l)p - 2b 
(n + l)ax'^ + bx + c 



dx 



and it can be assumed up to G^^i -equivalence (see Theorem II ip that, if n 7^ —1, the parameter tuple (a, b, c, d,p) 
takes only the following inequivalent values: 

{(0,l,0,d,(g + 2)/(n + m-l)), (0,0,l,l,p), (0,0,1,0,1), (l/(n + 1), 0, 1, 

where p is an arbitrary constant; d ^ and, if d = 0, p ^ 0; 

3n + 4 



{d,q) / (0,0), 



-,-3 



+ 4 , ..T 3n + 4 3 
d —— — and, it d = —— — — rv, Q ^ ^7: ^ 



n + 1' ' n + lj' 2(n+l) ' 2(n + 1) ' ^ " 2 2(n + 1) 

If n = —1, up to G^=i -equivalence the parameter tuple (a, 6, c, d,p) can be assumed to belong to the set 

{(0,l,0,d',p'), (0,0,1,0,1), (£,0,l,0,p")}, 

where d' and p' are arbitrary constants, p" ^ 0. In Case 5 the parameter p should be neglected. In Case 8 the 
parameter-functions / and h can be additionally gauged with equivalence transformations from G'^=i^rn=n+i (see 
Theorem I12|l . For example, we can put f — 1 if n ^ —4/3 and f — otherwise. 



there is a contraction of Case 9 of Table 2 to Case 8 of Table 1. The corresponding equations 
from classes ([3]) and ([2]) are 

(e™^is)- + ax-V", (35) 
= (n"n^)^ + Qx- (36) 

whose maximal Lie invariance algebras g and g are generated by the vector fields Xi = 5j, 
X2 = nidi ~ ^u-, X^ = nxdx + 2du and Xi = dt, X2 = ntdt — udu, X^^ = nxdx + 2udu-, 
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respectively. The contraction 2.9 — t- 1.8 can be realized using the simpler transformation 

t = t, X = X, u = 5{u — l), n=— , a = 5a (37) 



than transformation (j33p . In the course of this contraction the algebra q is contracted to the 
algebra g as a Lie algebra of vector fields in the space of {t,x,u). Namely, Xi — t- Xi, X2 — s- X2 
and X3 — )• X3. 

Let us study the related contractions of Lie reductions of equation (j36p to ones of equa- 
tion (|35p . Inequivalent Lie reductions of these equations with respect to one-dimensional subal- 
gebras of the corresponding maximal Lie invariance algebras are exhausted by those presented in 
Tables 3 and 4. For convenience we omit tildes in Table 4. The transformations of the invariant 
independent and dependent variables, which are induced by transformation (]37p . take the form 
(p = 6{ip — 1) and w = w in all cases of Table 3. 



Table 3. Lie reductions for Case 9 of Table 2. 



N 


X 




u — 


Reduced ODE 


1 


X2 - 




\t\- — ^{uj) 


{(p"(foj)uj - iJ-ij^eipui + ^^^S'fi + Qa;~^(^"+^ = 0, £ = sign* 


2 


Xi 


t 


\x\iLp{uj) 


riV^ - {an^ +2n + 4)v3"+^ = 


3 










4 




X 


(p{uj) 


/ n \ 1 —2 n+l r\ 

{if (ftojio -\~ aoj ip ^ =0 



Table 4. Lie reductions for Case 8 of Table 1. 



N 


X 




u = 


Reduced ODE 


1 


X2 - ^iX3 




^H-i±2iiln|t| 


(e"'^)„„ - Ainea;vp„ + (1 + 2^)^ + anoj-^e""^ 


= 0, e = signt 


2 


X-i 




V3(a;) + |ln|a;| 


nv9^ - (an + 2)e"*' = 




3 






± |t 


(e"'^)„„ ± na;(^^ =F 2 + antj-^e"'^ 


= 


4 




X 


V5(a;) 







Consider Case 1 of Table 3 in detail. The transformed version and the corresponding limit 
of the ansatz n = |t| ^ ^{oj) are 

Therefore, the contracted ansatz is = (p{(jj) — ^"^^^ In The reduced equation from Case 1 of 
Table 3 is mapped by transformation ([37P to the equation 



5h 



5n + l 



^-1 



99 \ an ( (p\ 



fifieuxfcj + (1 + 2/i)e (1 + t)+^(1 + t) =0- 



Then the limit process at 5 — )• -|-oo leads to the equation 

~ - an ~ ~ 

(e«^) - - - fineO^^ + (1 + 2f^)e + ^e"^ = 

which is also obtained from equation (|35p by the reduction with respect to the contracted ansatz 
and presented by Case 1 of Table 4. 
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Analogously we obtain contractions of reductions 3.2 — t- 4.2, 3.3 — t- 4.3 and 3.4 — )• 4.4. 

For Cases 2 and 4 of Table 3 exact solutions of reduced equations were found in [23]. The 
substitution of these solutions to the respective ansatzes results in the following exact solutions 
of equation (p6|) : 



X 



C - (an + 2 + 4n-i)t 



|CiV^lnx + C2^/i|"+i, 

1 



if a' = 0, 

Cix^i + C2X^2 1 — ^ if a' > 0, 

1 

Ciy^sin((7lnx) + C2-v/xcos((7lnx) , if a' < 0, 



(38) 



(39) 



where 



a 



l-4a(n + l), xi^2 



1 ± Va' 



a 



Here and in what follows C, C\ and C2 are arbitrary constants. Applying transformation (|37|) 
to solution p8]) and proceeding with the limit b — )• +cxd, we obtain 



1 + 7 



<5n 



^2 C 



4 \~ 

an + 2 + — k 
no/ 



[C - (dn + 2)ty 



As a result, we construct the exact solution 



M = — In 

h 



£2 



C - (on + 2) t 



for equation (j35]) . Applying the same technique to solutions (|39ll leads to the steady-state exact 
solutions of ([5^ : 



''i In |Ci\/ilnx + C2\/^|, 
ilnjCii^i + C72X^2|^ 
In sin((T Inx) + C2\/icos((T lnx)| 



if dl = 0, 
if a' > 0, 
if d' < 0, 



where 



a 



1 — 4an, xi 2 



1 ± Vd' 



Another way for finding this solution is to integrate the reduced equation of Case 4 from Table 4. 
By the obvious transformation = e"''' the reduced equation is mapped to the Euler equation 



6.3 Contractions of conservation laws 

We use contractions in order to construct conservation laws of equations from class ([3]) with g = 1 
using results obtained in [23] for equations from class ([2]) with the same gauge of g. Note that 
the consideration can be easily extended to the entire classes ([3]) and ([2]) using transformations 
from the corresponding equivalence groups. 

Roughly speaking, a conservation law of a system C of differential equations is a diver- 
gence expression that vanishes on solutions of this system. Thus, in the case of two indepen- 
dent variables t and x and one unknown function u the general form of conservation laws is 
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DtF{t, x,U(^j.)) + DxG{t,x,U(^r)) = whenever ti is a solution of C. Here Dt and are the 
operators of total differentiation with respect to t and x, respectively, and denotes the set of 
all the derivatives of the functions u with respect to t and x of order not greater than r, including 
u as the derivative of the zero order. The components F and G of the conserved vector {F, G) 
are called the density and the flux of the conservation law. Two conserved vectors {F, G) and 
(-F', G') are equivalent if there exist such functions F, G and H oi t, x and derivatives of u that 
F and G vanish for all solutions of £ and F' = F + F + D^H , G' = G + G - DfH . A conserved 
vector is called trivial if it is equivalent to the zero conserved vector. 

It is found in [23] that there are only two subclasses of equations of form ([2]) which admit 
nontrivial conserved vectors. Thus, assuming the gauge g = 1, each equation from class ^ 
with m = n + 1, i.e., an equation of the form f{x)ut = {u^Ux)x + h{x)u^~^'^^ admits two 
linearly independent conservation laws with the following conserved vectors (F* , ) and the 
characteristics A*, i = 1,2: 

n^-1: (^7n, -(^Vu, + (/.i.^), A^ = (^^i = l,2; 



n 



-1: {xfu, —xu ^Ux + hiu), A^ = x; {fu, —u ^u^), A^ = 1. 



Here /3i and /32 are arbitrary constants. The functions 93* = f^{x), i = 1,2, form a fundamental 
set of solutions of the second-order linear ordinary differential equation ip^x + {n + l)hip = 0. 

The other subclass of equations admitting nontrivial conserved vectors is singled out from ^ 
under the gauge g' = 1 by the conditions m = 1 and h = af, where a is an arbitrary constant, 
i.e., it consists of equations of the form 

f{x)ut = KuxOx + af{x)u. (40) 

The corresponding conserved vectors and characteristics are 



—nt 

xe 



' ' / —at f„. „—at„,n„, \ \2 „—at. ^ ' 



n = —1: 



( xe-'^'fu, e-'-\-xu^Ux + ^) ) , A^ 
(e-"7M, -e-°*u"n^), A^ = e 
xe~"*/M, e~^^{—xu~^Ux ^\^u)\ }^ = xe~°^^ 
' e-"Vu, -e-'^^u-^Ux ), A^ = e""*. 



In order to contract equations from class ([2]) to equations from class ([3]), we should vary 
the arbitrary element n. This is why only the case of general n is appropriate for contractions. 
There are three different ways in order to realize contractions of conservation laws. Namely, we 
can contract characteristics of conservation laws, their conserved vectors or conservation laws 
as divergent expressions themselves. 

We illustrate these possibilities in detail using equations (j40p with n 7^ — 1 and their conser- 
vation laws associated the same characteristic A^ = xe""^*. The corresponding conserved vectors 
are presented in (jHJ. At first we apply equivalence transformation (f37|l to equation (HO]l and 
proceed to the limit b — >• +00. As a result, we obtain the class of equations (tildes are omitted) 

/(x)txt = (e™n,)x + a/(x), (42) 



i.e., equations from class (llip with m = and h = af. 

For the image A^ of the characteristic A^ = xe~"* with respect to transformation ()37p we 
have that A^ — )• x if 5 — )• +00. Now we are able to construct the corresponding conservation 
law of ()42p using the characteristic obtained as an integrating factor. After the multiplication 
by X the equation (j42p can be written in divergent form as 



A (xfu - ax ft) + Dx ( -xe""u^ + -e"" ) = 0. (43) 



n 



21 



Therefore, we construct conservation law (|43p of equation (j42p via carrying out a limiting pro- 
cess of characteristics. Another way is to directly deal with divergence expressions. Thus the 
conservation law 



of (jiOj) with n / — 1 is transformed by (j37j) to 

A- (i + ) ^ (-"^'^^ ^ 0" T - {I - - °- 

Multiplying the obtained expression by b and adding the term —xfb under D^, we proceed 
to the limit b — )• +00. (As the term —xfb depends only on x, it is negligible in view of the 
differentiation with respect to t.) Omitting tildes, we exactly obtain the conservation law (|43p . 

Conserved vectors can be contracted in the same way using the fact that we can add expression 
which depends on x only to the density component up to equivalence of conserved vectors. 

Contracting the conservation laws obtained in [23] jointly with the corresponding equations, 
we obtain the following assertion. 



Theorem 10. A complete list of equations from class (jlip possessing nontrivial conservation 
laws is exhausted by the following ones. 

1. m = n : 

{ip'fu, -^ie^^u^ + ^ifie^^), ^\ i = l,2. 

2. m = 0, h = af : 

{xf{u-at), -xe""tx, + ie™), x; {f{u-at), -e""n,), 1. 

Here the functions (p^ = <p^{x), i = 1,2, form a fundamental set of solutions of the second-order 
linear ordinary differential equation (p^x + nhip = and a is an arbitrary constant. 

Simultaneously with constraints on the arbitrary elements we also present conserved vectors 
and characteristics of the basis elements of the corresponding space of conservation laws. 



7 Conclusion 

In this paper we carry out the extended Lie symmetry analysis of equations from class , where 
A and B are exponential functions of u. In this sense the present paper continues the series of 
papers [22H25] . where equations ([TJ with power nonlinearities and semilinear equations of this 
form with exponential source were studied. 

Lie symmetries, admissible transformations and conservation laws of equations from class ([3]) 
are exhaustively classified. This difficult task is achieved due to the usage of generalized equiv- 
alence groups instead of usual ones, nontrivial gauging of arbitrary elements and the repre- 
sentation of class (l3|) as a union of normalized subclasses. Moreover, in the course of group 
classification of equations with m = n the associated conditional equivalence group is used, 
which essentially simplified calculations and the final result. 

We also study limit processes, called contractions, between equations from classes ([2]) and ([3]) 
jointly with limit processes between the corresponding Lie invariance algebras. This allows us to 
predict and then to check the results of group classification for class ([3]). Moreover, we construct 
conservation laws and exact solutions for equations from class ^ using contractions and the 
related results obtained for equations of the form ([2]) in [23j. 
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The feature of this paper in comparison with the preceding ones is that the method of furcate 
spht apphed to classify Lie symmetries classification as well as the usage of the direct method 
for finding equivalence and admissible transformations are described in detail. 

Existing examples on group classification of variable-coefficient diffusion-convection equa- 
tions [Hl Um fTU] show that the subclasses with power or exponential nonlinearities A and B are 
usually the most complicated to be classified. In a forthcoming paper we intend to complete Lie 
symmetry analysis of equations from class ([T]), where the nonlinearities A or B are not power 
or exponential functions. 



A Gauging of parameters in classification list 
with power nonlinearities 

We briefly describe results on the group classification of class ([2]) obtained in [23], which are 
needed for the proper understanding of the final list of inequivalent Lie symmetry extensions 
collected in Table 2. In order to attain to the complete similarity to the presentation of the 
group classification of class ([3]) in the present paper, we carry out certain modifications and 
enhancements of results from [23]. In particular, we thoroughly gauge parameters appearing in 
Cases 2 and 5 of Table 2 and partition these cases into simpler inequivalent subcases in the same 
way as this is done for class ([3]) in Section [H 

Via gauging of the arbitrary element g by transformations from the corresponding equivalence 
group, the group classification of class ([2]) is reduced to that of its subclass singled out by the 
constraint g = 1. 

Theorem 11. The generalized equivalence group G'^^i of the subclass of equations having the 
form ([2]) with g = 1 consists of the transformations 

r ^ , ^ ~ + 67 ^ ^ , ^ ^ , i_ 

t = dit + d2, x = - —, n = ds c)4X + 

~ (5l(53", ,3n + 4 ~ ^3"-'"+^ , m+3n + 3 

f = ^p-\hx + 5^} f , h = — ^2 — \5iX + 55\ "+i h, n = n, m = m, 

where 5j, j = 1, . . . ,7, are arbitrary constants such that did^ ^ 0, A = S^Sq — 5461 7^ and the 
tuple {S4,,S5,Sq,6y) is defined up to a nonzero multiplier, e.g., we can set A = ±1. The arbitrary 
element n is assumed to be unequal to —1. For n = —1 transformations from the group G'^^i 
take the form 

i=6it + 62, x = 64X + 65, ii = 6se^^^u, 

f= T— ^-T — f: h = ; — /l, TL = U, 771 = 771, 

where 5j, j = 1, . . . ,6, are arbitrary constants, 7^ 0. 

Theorem 12. The conditional equivalence group C^^^^^^^^^i of class ([2|) associated with the 
constraints g = 1 and m = n + 1 is formed by the transformations 

t = 6it + 62, x = ip{x), u = ip{x)u, 

Here (p and ip are arbitrary smooth functions of x and 6j, j = 0,1,2, are arbitrary constants 
satisfying the conditions 5qlPx = ijP'^^'^ and 6o6iip / 0. 
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In the course of group classification in the general case m ^ l,n + 1 we derive that the 
maximal Lie invariance algebra of an equation from class ^ with (7 = 1 is a proper extension of 
the kernel algebra {dt) only if the corresponding value of the arbitrary elements / and h satisfy 
a system of the form 

fx 

((n + l)ax^ + bx + c)-^ = — (3n + 4)ax + d, 

L ^^^^ 

{{n + l)ax^ + bx -\- c)-^ = — (3(n + 1) + m)ax + (1 + n — m)p — 2b, 

where a, b, c, d and p are constants which are not simultaneously equal to zero. 

Lemma 2. If n ^ —1, up to C^^^- equivalence the parameter tuple {a,b,c,d,p) can be assumed 
to belong to the set 

{{0,l,0,d,p), (0,0,1,1,P), (0,0,1,0,1), (0,0,1,0,0), {l/{n + l),0,l,d,p)}, 
where d^ and, if d = 0, p ^ 0; p is an arbitrary constant; 

Ifn = —1, C^^i-inequivalent values of the parameter tuple (a, b, c, d,p) are exhausted by elements 
of the set 

{(0,l,0,d',p'), (0,0,1,0,1), (0,0,1,0,0), (e",0,l,0,p")}, 
where d' and p' are arbitrary constants, e" = ±1 and p" ^ 0. 

Proof. Combined with the multiplication by a nonzero constant, each transformation from the 
equivalence group is extended to the coefficient tuple of the system (I44p . The extension 

takes the form 

(n + l)d = v{5l{n + l)a - d^d^b + 5|c), 

b = u{-265S7{n + l)a + {6467 + dsde)^ - 26466c), 

3n + 4 

c = v{6j{n + l)a — 666ib + (5|c), d = uAd H --^i'{6567{n + l)a — 64,67b + 646QC), 

p = vAp z^((55(57(n + l)a — 6467b + 646QC) 

n + 1 

if n 7^ — 1 and 

d = —{a + 6Qb), b = vb, c = v{64C — 5^b), 

V (46) 
d= iyd+ —{65a + ^5^6^ - Sg^ac), p = up - —{65a + 6566b - 6664c) 
04 O4 

if n = —1. Here A = 6^66 — 6467 7^ and is an arbitrary nonzero constant. 

If n 7^ — 1, the proof is similar to the proof of Lemma [1] In this case there are only three 
Gj'-inequivalent values of the triple (a, b, c) depending on the sign oi D = b^ — Anac, 

(0,1,0) if D>0, (0,0,1) if D = 0, (l/(n + 1), 0, 1) if D<0. 

Indeed, if D > then there exist two linearly independent pairs {64,65) and {66,67) such that 
(5|(n + l)a — 6465b + 6lc = and 6j{n + l)a — 6667b + 6'^c = 0. For these values of (5's we have d = 
c = 0. In the case D = we choose values of 64, 65, 66 and 67 for which 5|(n+l)a— (54(556+(5|c = 
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and the pair {6q,6j) is not proportional to the pair {84^,55). Then we obtain that a = and 
b = i'5j{5/ib — 255(77, + l)a) + iy5Q{55b — 2840) = 0. The residual coefficient (6 if > and c if 
-D = 0) is necessarily nonzero and hence it can be scaled to 1 using the multiplication by the 
appropriate value of v. If D < 0, we have ac 7^ and can set a > 0. As the matrix 

/ (n + l)a -6/2 \ 

is symmetric and positive, the corresponding bilinear form is a well-defined scalar product. 
Choosing v = 1 and pairs [84, 5^) and {5q,5j) which are orthonormal with respect to this 
product, we obtain (n + l)a = c = 1 and 6 = 0. 

Now for each of the above inequivalent form for both the tuples (a, 6, c) and (a, 6, c) we look 
for possible gauges of the coefficients d and p. 

Thus, from {a,b,c) = (a, 6, c) = (0,1,0) we derive 6465 = 6^67 = and 6467 + 656^ = v^^ . 
This system in (5's has two solutions, 84 = 87 = Q with A = 8^8^ = and 8^ = 8q = Q with 
A = —8487 = The first solution leads to the identical transformation of the coefficients d 

and p. For the second solution the transformation takes the form d = —d — (3n + 4)/(n + 1), 
p = —p + l/(n + 1) . This is why up to G~^;^-equivalence we can assume that d and p satisfy (05]) . 

Setting (a, 6, c) = (a, 6, c) = (0, 0, 1) results in 84 = 0, v8l = 1 and A = 8^8^ 7^ 0. The 
transformation of d and p is reduced to simultaneous scaling with the same multipliers, d = z^Ad 
and p = v/\p. This allows us either to set d = 1 if (i 7^ or to scale p ii d = Q and hence d = 0. 
As a result, we obtain the tuples (0,0, 1,1, if d 7^ 0, (0,0,1,0,1) if d = and p 7^ 0, and 
(0, 0, 1, 0, 0) in the case d = p = Q. The last tuple corresponds to equations of the form ([2]) with 
constant coefficients (Case 3 of Table 2). 

The equality ((ra + l)a, 6, c) = ((n + l)a, 6, c) = (1, 0, 1) implies 81 + 8l = 8l + 8"^ = and 
848^ + 8587 = 0. Hence 8q = e^s and 87 = —684, where e = ±1. The transformation of the 
coefficients d and p is reduced to the multiplication hy e, d = ed and p = ep. This is why we 
can only set (i ^ and, if d = 0, p ^ 0. 

Consider the case n = — 1. As the tuple {a,c,b,d,p) is nonzero, system (06]) implies that 
(6, c) 7^ (0, 0). If the tuple (a, 6, c) satisfies the condition 6 7^ (resp. 6 = a = 0, resp. 6 = and 
a 7^ 0) then it is G~^]^-equivalent to the tuple (0, 1,0) (resp. (0,0, 1), resp. (e",0, 1)). Similarly 
to the case n 7^ — 1, now we look for possibilities of gauging the parameters d and p after fixing 
one of the above inequivalent forms for both the tuples (a, 6, c) and (a, 6, c). 

The equality (o, 6, c) = (a, 6, c) = (0,1,0) implies 1^ = 1 and 85 = 8^ = 0. Therefore, the 
parameters d and p are identically transformed. 

Setting (a, 6, c) = (a, 6, c) = (0, 0, 1) in (j46]l results in 84 = v^^, d = v{d—8Q) and p = ^{p+Sq). 
Therefore, we can set {d,p) to be equal (0,0) or (0, 1). 

It follows from (Ii6]) with (a,6, c) = (a, 6, c) = (e",0, 1), where e" = ±1, that u = 84 = ±1, 
d = vd + e"8^ — v8q and p = vp — e"8^ + v8q. This allows for setting one of the parameters d and 
p to zero. After fixing the zero value of d, we can additionally alternate the sign of p. Hence we 
assume {d,p) = {0,p"), where p" ^0. □ 

Lemma [2] implies that up to G~^]^-equivalence Case 2 of Table 2 is partitioned into three 
inequivalent subcases, 

1. {f,h) = {\x\'^,e\x\'^), q = {I + n — m)p — 2: {dt, {d + 2 — pn)tdt + xdx + pudu); 

2. {f,h) = {e'^^,ee'^^), q = {1 + n — m)p: {dt, {d — pn)tdt + dx + pudu); 

3. n7^ -1, (/,/l) = (^(x2 + l)-4STTyerfarctanx^g^^2^;L^-2f(t^gqarctanx^ ^ {l+n-m)p: 

{dt, {d - pn)tdt + (x2 + l)dx + (x/(n + 1) + p)udu)] 

n = —l, {f,h) = (^^'^"^^ ,ee^^"^^^'^^^ , q = —mp: {dt , qtdt — mdx + (q + me" x)udu) ■ 
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Here e,e" = ±1, {d,q) / (0,0) for the first and second subcases and additionally {d,q) / 
(— (3n + 4)/(n + 1), —3 — m/{n + 1)) in the first subcase with n ^ —1 and {d, q) ^ (0, 0) in the 
third subcase with n = —4/3; otherwise we have cases of Lie symmetry extensions of greater 
dimensions. It follows from Lemma [2] that up to -equivalence we can set certain constraints 
for the parameters d and q. (It is convenient to use q instead of p as a parameter.) These 
constraints are different for the cases n ^ —1 and n = — 1. 

If n 7^ —1, for the first subcase an exhaustive gauge implied by G~^^-equivalence consists of 
the inequalities 

3n + 4 3n + 4 ^ 3n + 3 + m 

a ^ r and, ii a = r, q ^ : — . 

2(n + l) ' 2(n + l)' ^ 2(n + 1) 

They can be set using the equivalence transformation 
1 , , i_ 

t = t, X = —, -U = X 
X 

whose extension to the parameters d and q is given by d = —d — (3n + 4)/(n + 1) and q = 
—q — 3 — m/(n + 1). In the second subcase the parameters d and q can be gauged using a scaling 
of X. More precisely, d = 1 mod G'^^^ if d 7^ and q = 1 mod G'^^-^ if d = 0. In the last 
subcase we can just simultaneously alternate the signs of d and q. Hence the exhaustive gauge 
is presented by d ^ and, if d = 0, g ^ 0. 

Consider n = — 1. Then in the first subcase the parameters d and q are not changed by 
transformations from G^^i- In the second subcase the parameters d and q can be gauged and 1, 
respectively. In the last subcase we can alternate the sign of q and assume g ^ mod G~^]^. 

A similar partition can be also carried out for Case 5 of Table 2. Namely, we have the 
following inequivalent subcases: 

1. f = \x\'^: {due-'''\dt + eudu),nxd^ + {d + 2)udu)\ 

2. f = : (at, e-^"*(at + eudu), nd^ + ud^)- 

3. n / -1, / = {x^ + l)-2fST) grfarctanx . 

{du e-^''\dt + eudu), n{x'^ + 1)9^ + (nxjin + 1) + d)udu); 

n = -l,f = e¥'-' : {du e'"\dt + e"ud^), 3^ + (e" x - d)ud^). 

In the first subcase d ^ and, if n 7^ —1, d ^ — (3n + 4)/(n + 1). In the third subcase d ^ 
if n = —4/3 and e" = ±1 if n = —1. The gauges for the parameter d modulo G^^^^-equivalence 
coincide with the gauges of d in the respective subcases of Case 2 from Table 2. 
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